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AN ANALYTICAL STUDY OF TURBULENCE RESPONSES, 

INCLUDING HORIZONTAL- TAIL LOADS, OF 
A CONTROL-CONFIGURED VEHICLE WITH 
RELAXED STATIC STABILITY* 

Boyd Perry HI 
Langley Research Center 

SUMMARY 

A static stability augmentation system (SSAS) with angle-of-attack, pitch-rate, and 
forward-speed feedback is used to stabilize control-configured rigid vehicles with relaxed 
static stability. Airplane configurations include a baseline (statically stable) configura- 
tion and two representative control- configured vehicle (statically unstable) configurations. 
The first statically unstable configuration has baseline geometry and an aft center of 
gravity; the second statically unstable configuration has reduced tail length and area and 
a forward neutral point. Stability, flying-qualities, and maneuverability requirements are 
imposed on the airplane-SSAS systems. Turbulence responses of the rigid airframe and 
tail loads are examined in terms of configuration changes for a given control law and in 
terms of control law changes for a given configuration. Results indicate that even though 
the tail loads for the unstable configuration with the small tail are less than those for the 
stable configuration, the tail-load intensity is larger. 

INTRODUCTION 

Analytical studies of turbulence responses of autopilot- controlled vehicles that have 
been reported in the open literature include parametric variations in feedback gains 
(refs. 1 and 2), airplane geometry and flight conditions (ref. 2), and Mach number (ref. 3). 
In each study, the autopilots maintained the example airplanes on prescribed flight paths 
and the example airplanes were inherently longitudinally statically stable at the flight con- 
ditions investigated. In an earlier study (ref. 4) the elastic- and rigid-body turbulence 
responses of a long-range bomber airplane were calculated for a number of stability aug- 
mentation systems. Again, the airplane was inherently longitudinally statically stable. 

*Part of the information presented in this report was included in a thesis entitled 
"An Analytical Study of Turbulence Responses, Including Horizontal-Tail Loads, of a 
Control-Configured Jet Transport With Relaxed Static Stability" submitted to The George 
Washington University in partial fulfillment of the requirements for the degree of Master 
of Science, February 1975. 



The recent emergence of the control -configured vehicle (CCV) concept has introduced air- 
planes which are deliberately designed to be loi^tudinally statically unstable (referred to 
as relaxed static stability) at certain flight conditions within their flight envelopes. (See 
refs. 5 to 10.) These vehicles rely on static stability augmentation systems (SSAS) for 
longitudinal static stability augmentation at all times when the vehicle is in flight. (See 
refs. 9 to 11.) Because the CCV concept is new, there are some technical areas which 
have received little attention in the open literature. One such area is the study of control- 
configured vehicles with relaxed static stability flying through continuous atmospheric 
turbulence. 

The purpose of this paper is to determine airframe and horizontal-tail-load turbu- 
lence responses of a control -configured rigid vehicle with relaxed static stability (RSS). 
The vehicle is representative of small executive jet transports. Three variations of the 
same airplane configuration are investigated at two different flight conditions. The three 
airplane configurations include a baseline (longitudinally statically stable) conventional 
configuration and two configurations with relaxed static stability (one with an altered cen- 
ter of gravity and one with reduced tail length and area). The configurations with RSS are 
artificially stabilized with a static stability augmentation system which features angle-of- 
attack feedback, pitch-rate feedback, and forward-speed feedback to the elevator. The 
flight conditions are Mach numbers of 0.50 and 0.75 at an altitude of 6100 meters. Random 
process theory is used to determine the rigid-airframe and tail-load turbulence responses. 
Only the vertical component of turbulence is considered. The turbulence varies in the di- 
rection of flight but does not vary along the span (one-dimensional gust field) and it is rep- 
resented by the Von Karman power spectrum. The rigid-airframe responses include the 
longitudinal degrees of freedom and other derived quantities; and the tail-load responses 
include torque and bending moment of the horizontal tail. Stability requirements, flying- 
qualities requirements, and maneuverability requirements are imposed on the SSAS in each 
configuration with RSS. These requirements are first satisfied by appropriate feedback 
gain selections and then the turbulence responses of each configuration with RSS are calcu- 
lated and compared with the turbulence responses of the baseline configuration. In addi- 
tion to the analysis and discussion of the turbulence responses given in the paper, two 
appendixes are presented: Appendix A defines flying-qualities criteria in the complex 
plane; appendix B defines the apparent maneuverability for the augmented airplane and 
examines the effect of the feedback gains on the apparent maneuverability. 

SYMBOLS 

A(s) aerodynamic coefficient matrix 

Aj^ ratio of root- mean- square value of response R to root- mean- square value 

of gust input 


2 



b 

bt 

Cl 

Cl,o 

CM,b 

Cm 

^mq 

Cma 


Ct 


Cx 





c; 


a 


c 


Ct 


Ct 


f(s) 

g 


wing span, m 
horizontal-tail span, m 
tail-load coefficients (see eqs. (8) and (9)) 
lift coefficient 

lift coefficient for steady flight condition 
bending- moment coefficient, 

pitching-moment coefficient 
apparent stability derivative 

apparent Cm^^ stability derivative 

T 

torque coefficient, = 

x-force coefficient 
z-force coefficient 
apparent stability derivative 

apparent Cz^y stability derivative 
wing mean aerodynamic chord, m 
local chord of horizontal tail, m 
horizontal-tail mean aerodynamic chord, m 
vertical-gust column matrix 
acceleration due to gravity 


L 



frequency-response function of response R 
h center-of-gravity position in terms of c 

hjjj stick-fixed maneuver point in terms of c 

apparent stick-fixed maneuver point in terms of c 
hj^ stick- fixed neutral point in terms of c 

hn apparent stick-fixed neutral point in terms of c 

lyy airplane pitch moment of inertia, l^-m2 

ig nondimensional pitch moment of inertia, Slyy^^Sc^ 

Ku forward-speed feedback gain, rad/m/sec 

Kq, angle-of-attack feedback gain, rad/rad 

Kg gear ratio, rad/rad 

pitch-rate feedback gain, rad/rad/sec 
kj^ undamped natural reduced frequency, (v^cj2\XQ 

kQ damped natural reduced frequency, ct)QCy^2uo 

L scale of turbulence, m 

distance from center of gravity to horizontal-tail aerodynamic center, m 
M moment about Y-axis; also Mach number 

horizontal-tail bending moment outboard of 24 percent semispan, N-m 
m mass of airplane, kg 

n acceleration factor, g units 
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vertical acceleration of airplane center of gravity in g units, positive down 

column matrix of independent variables 

general turbulence response (u, An, M'd, etc.) 

wing area, m2 

horizontal-tail area, m2 

Laplace variable 

horizontal- tail torque about 40 percent chord, N-m 

airplane characteristic time, c/2uq 
elevator servo characteristic time, sec 
forward-speed perturbation, m/sec 
nondimensional forward-speed perturbation, u/uq 
airplane forward speed, m/sec 

perturbation velocity in z-direction, positive down, m/sec 

vertical-gust velocity, positive up, m/sec 
force in x-direction 
force in z-direction 

airplane angle-of-attack perturbation, w/uq 
gust angle of attack, -w^^Uq 
horizontal -tail angle of attack 

elevator deflection perturbation, positive trailing edge down, rad 

elevator deflection necessary to hold a pull-up maneuver of load factor n, rad 
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A5’ 
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C 

d 

At 

A.25 

M 

P 


a 





T 


^Wg(‘^) 

OJ 




n 


Ct) 
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apparent elevator deflection necessary to initiate a pull-up maneuver of load 
factor n, rad 

control-stick deflection, rad 
downwash angle at horizontal tail, rad 
damping ratio 

pitch-angle perturbation, positive nose up 

sweep angle of horizontal-tail 40 percent chord, deg 

sweep angle of wing quarter chord, deg 

nondimensional airplane mass, 2m/pSc 

air density, kg/m^ 

real part of complex number 

root- mean- square value of response R 

root- mean- square value of vertical-gust velocity, m/sec 

transport time lag, S-t/nQ 

power spectral density function of vertical-gust velocity, (m/sec) 2yrad/sec 

circular frequency, rad/sec 

undamped natural circular frequency, rad/sec 

damped natural circular frequency, rad/sec 
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Subscripts: 


Ciuq! " 


9Cm 




8Cm 



_ 8Cna 



9Cl 

8Q! 


Mjj bending moment 

sp short period 

T torque 

t horizontal tail 

u forward speed 

u forward speed (nondimensional) 

a angle of attack 

0 !|. horizontal -tail ar^le of attack 

An center-of-gravity vertical acceleration 
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elevator deflection 


9 pitch angle 

9 pitch rate 

Abbreviations : 

CCV control- configured vehicle 

FQC flying -qualities criteria 

RSS relaxed static stability 

rms root mean square 

SSAS static stability augmentation system 


A dot over a quantity represents the derivative of that quantity with respect to time. 

ANALYSIS 


Equations of Motion 

The equations of motion used in this study are the three degree-of-freedom, longi- 
tudinal, small-perturbation equations (ref. 12), and an equation which describes the 
behavior of the static stability augmentation system. The equations of motion are written 
in the stability axes system and are presented in equation (1). 


2Mt*s - Cx 


2Ct - C, 
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2/it*s - C 2 




1 - e 
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J|-2m - Czqjf 
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|l-e~ 


KuUq 

1 + tjS 


1 + tjS 


.Bt*2s2 


■Cmqt*s 


KgS 
1 + tjS 


-Cz, 


^Zq! 


1 - 


TS 



- Cmqjt*5 


1 - e 


Kg 

1 + tjS 


(1) 
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The stability derivatives are defined in the list of symbols; t* is the airplane character- 
istic time defined in the list of symbols; tj is the servo -actuator characteristic time 
discussed below; t is the transport time discussed below; fi and ig are nondimen- 
sional airplane mass and pitch moment of inertia, respectively; Kq, K^, and are 
feedback gains; u, a, 9, and 6 are the small-perturbation independent variables; 
ttg is the gust angle of attack; and s is the Laplace variable. 

Equation (1) contains terms which account for the transport lags in downwash from 
wing to tail and the transport lag in gust encounter from wing to tail. These terms are 
indicated by the rectangular boxes in equation (1) and they reduce to unity if the first two 
terms of power series representation are substituted for the exponential (ref. 12). The 
time T in these terms is the time required for the airplane to traverse a distance equal 
to that between its center of gravity and the aerodynamic center of its horizontal tail. 

The power series representation was used to calculate the roots of the characteristic 
equation (for example, in root locus studies). The exponential was used to calculate the 
turbulence responses from the power spectral density functions. 

The last row of equation (1) describes the elevator deflection for the static stability 
augmentation system and may be expressed as 


6 = 


1 

1 + tjs 




( 2 ) 


Equation (2) states that the elevator is commanded to deflect in proportion to u, a, oig, 
and 6 and it is assumed that sensors are available to measure u, a + 0 !g, and 9 
directly. It is also assumed that the a-sensor is located at the airplane center of gravity, 
pother locations will produce changes in the stability characteristics and in the turbulence 
responses. For example, if the a-sensor is located at the nose of the airplane, a degree 
of gust anticipation is introduced relative to the present system. This anticipation may 
compensate, to a certain extent, for the lag in the servo actuator. The nose location will 
also cause the a-sensor to respond to pitching; thus, changes in will produce 

changes in an effective Ka. ) The quantity - — in equation (2) represents the 

^ / 1 + tjS 

dynamics of the servo actuator and has one-minus-exponential type response to a step 
input. For such a response, tj^ is the time it takes the elevator to sweep through 
approximately two-thirds of its final deflection. In this study, tj = 0.066 sec (less than 
1/10 sec as recommended in ref. 13), and consequently very much less than the period of 
the short-period mode. Equation (2) represents an idealized situation in that elevator 
aerodynamic hinge moments do not affect the response of the elevator. 
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The airframe frequency-response functions describe the responses of the airplane 
to a unit sinusoidal one-dimensional gust field of varying frequencies. Equation (1) may 
be written as 



(3) 


The frequency-response functions are obtained by solving equation (3) for the independent 
variables and substituting iu) for s. When -Wg^UQ is substituted for a^, the solu- 
tion of equation (3) is 


where 




[h - (ia,)] 



(io>) 

< > 

H Q (iu>) 

H 5 (io)) 



(4) 


(5) 


Equation (5) is the matrix of frequency-response functions in per-unit-gust-velocity 
form. The frequency-response functions of center -of-gravity vertical acceleration and 
horizontal -tail total angle of attack are expressed in terms of the frequency -response 
functions of the independent variables as 




IW 


Uo| 


H ^ (io)) - H ^ (iw) 


Vic 


w 


g 


( 6 ) 




1 - e 


■iioT 8e_ 
90 !, 


|H^ (iw) + iw ^ H Q (iw) 


Uo 


9e \ -iwT 

da) 


(7) 
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Tall- Load Equations 

The tail-load geometry is illustrated in sketch (a). 



The horizontal tail is swept and tapered with a full-span elevator and a frictionless hinge. 
The quantities T and Mj^ are the torque and bending moment calculated about the two 
perpendicular axes indicated by the dashed lines. The torque axis is the 40-percent chord 
line. The bending-moment axis intersects the torque axis at the 24-percent semispan 
position. Torque and bending moment are calculated about their respective axes by sum- 
mation of the forces acting outboard of the bending- moment axis. 

Both torque and bending moment contain components due to aerodynamic and inertia 
forces. The aerodynamic components result from the motion of the vehicle, the atmo- 
spheric turbulence, and the deflection of the elevator and they are obtained by integrating 
the aerodynamic-force distributions over the tail surface. The inertia components restilt 
from the vertical plunging and pitching accelerations at the tail and they are obtained by 
integrating the mass distribution of the tail over the tail surface and multiplying the 
result by the acceleration at the tail. 
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The following tail-load equations are derived in a manner similar to that described 
in appendix B of reference 1; 


T = ^13 ^14^ (®) 

B = C 2 iQ!t + C 225 + C 23 An + C 246 (9) 


where is the horizontal-tail total angle of attack, 6 is the elevator deflection angle, 
An is the vertical acceleration at the airplane center of gravity in g units, and 6 is the 
pitching acceleration. The Qj coefficients are constants obtained from the integra- 
tion of the aerodynamic and mass distributions. Coefficients and Cj2 coeffi- 

cients of the aerodynamic terms, and coefficients Cj 3 and €^4 are coefficients of the 
inertia terms. 

Torque and bending moment are expressed as dimensionless tail-load coefficients 
by dividing equations ( 8 ) and (9) by the quantity ipug^Stc^. 


Ct - Y 


|puo2StCt 


Mv 


'M,b 


ipuo2stCt 


( 10 ) 


( 11 ) 


1 2 - 

Because the quantity -^pug S^c^- is proportional to the volume of the horizontal tail, 

Ci 

equations ( 10 ) and ( 11 ) are measures of the horizontal-tail loading intensity (tail loads 
per tail volume). Equations (10) and (11) are useful because they permit comparisons of 
tail loads on tails of different sizes. 

Frequency -response functions of the torque and bending moment of the horizontal 
tail are expressed in terms of the frequency -response functions of the independent vari- 
ables and equations ( 6 ) and (7) as 


H^(ico) = CiiH^^(ico) + Ci2H^(icc) + Ci3H^(io.) - Ci4w2Hj_(ia.) (12) 

^ ^ ^ Wg Wg 


^ C2iH^^(icj) + C22H 5 (ia>) + C 23 H^(io)) - C 24 Cq 2 h g (ico) 

^ ^ ^ ^ 


w„ 


( 13 ) 
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Frequency-response functions of the torque and bending- moment coefficients are 
expressed in terms of equations (12) and (13) as 


\ 

H -p (ico) 



(icj) = 




H 


Hp (i^*^) — — 

^M,b 1 


M,b 

Wg 


(ico) 


w, 


g 


fpUo^StCt 


( 14 ) 


(15) 


Statistical Representations 

The airplane is treated as a linear system and its responses to random atmospheric 
turbulence are determined in a statistical sense by use of random process theory. In 
random process theory the power spectral density function (power spectrum) contains all 
the required statistical information describing a Gaussian random process with zero 
mean. 


Atmosphe ric turbulence .- The Von Karman power spectrum (ref. 14) was chosen to 
describe the vertical component of one-dimensional atmospheric turbulence 


4>wcr(^) ^ 

TTUr 




nil/6 


1 + 1.339 ^ (jo 


(16) 


The quantity a; is the circular frequency in radians per second; Uq is the forward 
speed of the example airplane; L is the scale of turbulence, chosen to be 762 m; and 
^Wo- root-mean- square (rms) value of the vertical-gust velocity, assigned to be 

o 

1 m/sec. The Von Karman power spectrum approaches a constant value at the low cir- 
cular frequencies and is characterized by a -5/3 power law at the high circular 
frequencies. 

Airplane responses .- The root- mean- square values of the airplane responses are 
the quantities of interest in this paper. The rms value of a general response (normalized 
to the unit rms vertical -gust velocity) is calculated in the followii^ manner: 
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where Hj^(ia)) is the frequency-response function of a general response R, and (oj) 

is the Von Karman power spectrum of equation (16). In actual practice the upper limits 
of integration in equation (17) were between 100 and 200 rad/sec. As stated in the dis- 
cussion of the equations of motion, the full exponentials (within the boxes in eq. (1)) were 
retained in evaluating the integral in equation (17). Had the first two terms of the power 
series representation been used instead of the full exponential, the integral in equation (17) 
would not have converged. Equation (17), when used in conjunction with the frequency- 
response functions of equations (5) to (7) and (12) to (15), yields the rms values of the 
turbulence responses. For the remainder of this paper, the ratio will be 

replaced by the symbol Aj^. 

STATIC STABILITY AUGMENTATION SYSTEM 

The static stability augmentation system (SSAS) provides stability and adds damping 
to those configurations with relaxed static stability. The SSAS contains three feedback 
loops: an angle-of-attack loop, a pitch-rate loop, and a forward-speed loop, all of which 
act through the elevator. Each feedback loop has a corresponding feedback gain. Refer- 
ence 7 states that such feedback to a pitching- moment control surface is adequate to pro- 
vide stability to airplanes with relaxed static stability. A block diagram showing the air- 
frame dynamics, the SSAS feedback signals, and the atmospheric turbulence is shown in 
figure 1. 

SSAS Requirements 

Various requirements are imposed on the static stability augmentation system such 
as providing stability, satisfying flying-qualities criteria, and increasing maneuverability. 
These requirements are met by varying the SSAS feedback gains. 

Stability exists when all roots of a system characteristic equation lie in the left half 
of the complex plane. For the configurations with relaxed static stability, the originally 
statically unstable mode is stabilized by varying the SSAS gains until all roots are in the 
left half plane. The presence of the servo- dynamics transfer function in figure 1 results 
in an additional root in the system characteristic equation. This root (referred to as the 
servo root) is real and it must also lie in the left half plane for stability. 
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Flying qualities criteria (FQC) are those for level 1 as defined in reference 15 in 
terms of phugoid and short-period damping ratios and undamped natural frequencies. 
These FQC may be visualized as regions in the complex plane (referred to here as FQC 
regions) within which the phugoid and short-period roots lie (ref. 16). The FQC are sat- 
isfied by varying the SSAS gains until the phugoid and short-period roots lie within their 
respective FQC regions. Figure 2 shows the phugoid and short-period FQC regions in 
the complex plane. Appendix A contains the derivation of the form of the FQC regions as 
shown in figure 2. Because of symmetry about the real axis, only the portions of the FQC 
regions in the second quadrant are presented. The insets show the values of undamped 
natural frequencies and damping ratios which define each region. The fan shape of the 
short-period region indicates that short-period FQC are specified in terms of both fre- 
quency and damping. The shape of the phugoid region indicates that phugoid FQC are 
specified in terms of damping only. Because the FQC regions lie totally within the left 
half of the complex plane, stability exists whenever the FQC are satisfied. 

Defining flying-qualities criteria is complicated by the presence of the servo root. 
The phugoid and short-period flying qualities are determined by the locations of the phu- 
goid and short-period roots in the complex plane. If the servo root is not "close" to 
either the phugoid pair or the short-period pair of roots, it will not seriously affect the 
physical interpretation of those roots. However, if the servo root is close to either pair, 
it will cause coupled motion between a conventional mode and the servo mode; and the 
FQC are not defined for such motion. 

Maneuverability in the longitudinal degrees of freedom usually is defined by the 
elevator angle per g, which is the elevator deflection required to trim the airplane in 
steady accelerated flight (such as a steady pull-up maneuver) of constant load factor. 

This deflection is related to the control-stick deflection through the gear ratio. For con- 
figurations with relaxed static stability, the SSAS provides a portion of the required ele- 
vator deflection in steady accelerated flight. In such configurations an apparent maneu- 
verability is defined in terms of the control-stick deflection. The control-stick deflection 
corresponds to an apparent elevator deflection (through the gear ratio) which the pilot 
commands in conjunction with the SSAS to achieve the steady pull-up maneuver. This 
concept is mathematically described in terms of stability derivatives and SSAS feedback 
gains in appendix B. For the present study, the maneuverability requirement for config- 
urations with relaxed static stability is that the apparent elevator angle per g is 25 per- 
cent less than the baseline elevator angle per g. 

Control Laws 

The static stability augmentation system is configured to perform two alternate 
functions: (1) to match the flying qualities of the baseline short-period mode without 
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regard for changes in maneuverability, or (2) to satisfy the maneuverability requirement 
and still maintain satisfactory flying qualities as defined by figure 2. Certain combina- 
tions of feedback gains are chosen to perform each of these two functions and such com- 
binations of gains will be referred to as control laws. One control law to perform func- 
tion 1 and two control laws to perform function 2 are discussed below. 

Control law 1 uses and feedback gains to perform function 1 and unique 

combinations of these gains are found to exist. Because the short-period mode is insen- 
sitive to changes in the Ku feedback gain, Ky is set to zero. 

Control law 2 uses and feedback gains (Ky set to zero) to perform func- 

tion 2. Unique combinations of these gains do not exist but single pairs of gains are 
selected in the following manner: the smallest Kq, is chosen to satisfy the FQC and 
then is adjusted until the maneuverability condition is satisfied. 

Control law 3 uses Ky and feedback gains (K^ set to zero) to perform func- 

tion 2. Unique combinations of these gains do not exist either but single pairs of gains 
are selected in the following manner: is chosen to satisfy the maneuverability con- 

dition and then the smallest is chosen to satisfy the FQC. 

AIRPLANE CONFIGURATIONS AND FLIGHT CONDITIONS 

Three airplane configurations are examined in this study: one with inherent static 
stability and two with relaxed static stability (RSS). The configurations are intended to 
represent one conventional design and two possible CCV-type designs of the executive-jet- 
transport class of airplanes. The conventional, or baseline, configuration has nominal 
weight and geometry and a positive static margin. Both configurations with RSS, through 
either weight redistributions or geometry alterations, have negative static margins at the 
flight conditions investigated. The first of these has the same geometry and weight as 
the baseline configuration, but a negative static margin due to an aft center-of-gravity 
position. The second configuration with RSS has altered geometry and reduced weight: 
the tail area is reduced 50 percent and the tail length is reduced 10 percent; thus there is 
a 5-percent reduction in total airplane weight. In this case the negative static margin is 
due to a forward-neutral point. Figure 3 shows the relative sizes of the wing and tail and 
the relative locations of the center of gravity and neutral point for each configuration. 
Table I presents pertinent physical dimensions for each configuration. References 8, 10, 
and 11 indicate that the geometry alterations and weight reductions pertaining to the sec- 
ond configuration with RSS are possible for a control-configured transport aircraft. 

The tails of the configurations with RSS are assumed to have enough authority to 
satisfy control -moment requirements such as nose wheel unstick. Performing detailed 
tail-sizing and weight-distribution studies is beyond the scope of this paper, but such 
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studies must be performed when designii^ an actual aircraft. For simplicity, when 
referring to either of the two configurations with RSS, they will be identified by their 
geometry and static stability characteristics: nominal unstable for the first and reduced 
unstable for the second. 

The flight conditions investigated are Mach numbers of 0.50 and 0.75 at an altitude 
of 6100 m. Table n contains center-of-gravity position and stability parameters for all 
configurations at both flight conditions. Table HI contains the stability derivatives for 
all configurations at both flight conditions. 

RESULTS AND DISCUSSION 

This section of the paper contains four parts: the first three examine the dynamic 
characteristics (including stability, flying qualities, and maneuverability) and turbulence 
responses for each control law; and the fourth examines the turbulence responses for a 
given control law in terms of configuration changes and the turbulence responses for a 
given configuration in terms of control-law changes. 

Control Law 1 

Control law 1 uses angle-of-attack and pitch-rate feedback to match the flying qual- 
ities (that is, frequency and damping) of the baseline short-period mode without regard 
for changes in maneuverability. Angle-of-attack feedback provides stability to the RSS 
configurations and pitch-rate feedback provides additional damping. 

Dy namic characteristics .- Table IV contains the values of the and feed- 

back gains for control law 1, the values of the natural frequencies and damping ratios of 
the phugoid and short-period modes, and the values of the apparent elevator angles per g. 
Table IV confirms that the flying qualities of the baseline short-period mode are matched; 
at both Mach numbers the short-period natural frequencies and damping ratios of the RSS 
configurations match those of the baseline configuration almost identically. It also indi- 
cates that larger amounts of and are required for the reduced-unstable con- 

figuration than for the nominal-unstable configuration because the smaller tail of the 
reduced-unstable configuration results in less elevator effectiveness. Smaller amounts 
of and are required at the higher Mach number because aerodynamic forces 

are greater for a given elevator deflection. The last column in table IV contains the 
apparent elevator angle per g for each configuration and flight condition. The apparent 
elevator angle per g for the nominal-unstable configuration is within 2 percent of the 
value of the actual elevator angle per g for the baseline configuration. The apparent 
elevator angle per g for the reduced-unstable configuration is (within 4 percent) twice 
that of the baseline configuration, and reflects the 50-percent reduction in tail area. 
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Figure 4 shows the phugoid and short-period roots (corresponding to the phugoid 
and short-period parameters of table IV) plotted in the complex plane. The cluster of 
short-period roots for each Mach number indicates that the frequency and damping of the 
baseline short-period mode have been matched. Both the phv^oid and short-period roots 
are within their respective FQC regions. 

Figure 5 contains a typical root locus of the characteristic equation for changes in 

and K^. The reduced-unstable configuration at Mach number 0.50 was chosen for 
the example; however, the other configurations show the same trends. The solid lines 
represent increases in ^with = o) and the dashed lines represent increases in 

Kg ^with = O.SS). Figure 5(b) is an enlarged view of the area near the origin in 
figure 5(a). When both gains are zero, the airplane is statically unstable and the roots of 
the characteristic equation are indicated by x. As increases, the dead-beat (non- 

oscillatory) unstable root crosses into the left half plane and combines with one of the 
roots from the classical third mode (ref. 14) to form the augmented phugoid mode. The 
other root from the third mode combines with another originally dead beat but stable root 
to form the augmented short-period mode. With set at 0.88, is increased until 

the frequency and damping of the augmented short-period mode match those of the base- 
line short-period mode. There is an additional dead-beat root associated with the eleva- 
tor servo transfer function. The root is stable and it is located far to the left beyond the 
part of the real axis shown in figure 5(a). Its characteristic motion is primarily in pitch 
with a time to half amplitude of 0.049 sec, only slightly slower than the zero-gain time to 
half amplitude of 0.046 sec. 

Turbulence responses .- Table V contains the turbulence responses (the A’s) for 
all configurations and flight conditions for control law 1, and the subscripts on the A’s 
indicate the appropriate response. When examining some of the turbulence responses of 
the baseline configuration at both Mach numbers, it was observed that: (1) the angle of 
attack Aa varied inversely with Mach number; (2) the vertical component of center-of- 
gravity acceleration A^n increases with Mach number, which reflects the correspond- 
ing increases in forward speed and airplane lift-curve slope; and (3) both horizontal- tail 
torque and bending moment, A>p and Aj^^b, increase with Mach number and reflect the 
increases in the aerodynamic and inertia components of both. 

An examination of the turbulence responses from configuration to configuration at a 
given Mach number revealed that the angle-of-attack response remains practically un- 
changed at each Mach number. Angle of attack is primarily a short-period phenomenon. 
Because short-period flying qualities are matched almost exactly and because a -feedback 
predominates over ^-feedback, A(^ remains unchanged. The pitch angle Ag is smaller 
for configurations with relaxed static stability for both Mach numbers (reductions ranging 
from 20 percent to 45 percent). However, center -of-gravity vertical acceleration A^j^ 
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increases between 10 percent and 15 percent in these configurations, probably as a result 
of the decrease in Ag. Elevator deflection Ag increases from the nominal-unstable 
to the reduced-unstable configuration because of the increased feedback gains for the 
reduced-unstable configuration. Horizontal -tail torque and bendir^ moment are larger 
for the nominal-unstable configuration than for the baseline configuration but smaller for 
the reduced-unstable configuration. The larger A values for the nominal -unstable con- 
figuration are primarily a result of the increase in the aerodynamic components from 
elevator deflection. Even though the tail loads are reduced (torque over 60 percent; 
bendii^ moment over 15 percent) in the reduced-unstable configuration, the horizontal- 
tail loading intensity (as indicated by the tail-load coefficients) is still larger at both 
Mach numbers. 

Inspection of the components of torque and bending moment revealed that approxi- 
mately 90 percent of each is due to aerodynamic contributions and that the remaining 
10 percent is due to the inertia contributions. Although the total aerodynamic forces on 
the tail due to and 6 are of the same order of magnitude, the net contribution to 
torque from elevator deflection is small compared with the net contribution from tail 
angle of attack. The aerodynamic force distribution due to elevator deflection has approx- 
imately equal contributions forward and aft of the 40-percent chord and results in consid- 
erable cancellation of torque. 


Control Law 2 

Control law 2 uses angle-of-attack and pitch-rate feedback to satisfy the maneu- 
verability requirement and at the same time maintain satisfactory flying qualities. 

Dy namic characteristics .- Table VI contains the values of and feedback 

gains for control law 2, the values of the natural frequencies and damping ratios of the 
phugoid and short-period modes, and the values of the apparent elevator angles per g. 
Table VI confirms that the maneuverability requirement (decreasing the apparent elevator 
angles per g of the RSS configurations 25 percent below the actual elevator angle per g 
of the baseline configuration) has been satisfied: at both Mach numbers the apparent ele- 
vator angles per g are within 1 percent of their required values. Table VI indicates 
that at both Mach numbers the short-period mode is overdamped (damping ratio greater 
than unity). It also indicates that larger amounts of are required for the reduced- 

unstable configuration because its smaller tail results in less elevator effectiveness. 
Smaller amounts of are required for the reduced-unstable configuration because the 
apparent elevator angle per g is a function of both and (see appendix B) and 

because the corresponding values of are larger. 

Figure 6 shows the phugoid and short-period roots (corresponding to the phugoid 
and short-period parameters of table VI) plotted in the complex plane. The phugoid and 
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short-period roots are within their respective FQC regions and thus confirm that the 
flying-qualities criteria have been satisfied. The overdamped short-period mode of the 
nominal -unstable configuration appears as two dead-beat roots at each Mach number, 
with corresponding motion primarily in plunge. Figure 7 contains a typical root locus 
for changes in and for the reduced-unstable configuration at a Mach number 

of 0.50. The present root locus initially proceeds in the same manner as the root locus 
for control law 1 in figure 5. However, in the present root locus, the vertical segment of 
the short-period mode stops just after entering the FQC region and, as a result, there is 
a lower damped natural frequency and a higher damping ratio. 

Turbulence respo nses.- Table VII contains the turbulence responses for all config- 
urations and flight conditions for control law 2. When examining the turbulence responses 
from configuration to configuration at a given Mach number, it was observed that the 
forward- speed response of the nominal-unstable configuration reflects the increased 
damping of the phugoid mode. For a Mach number of 0.50, Ajj for the nominal-unstable 
configuration is more than an order of magnitude less than for the baseline config- 
uration. For a Mach number of 0.75, A^ for the nominal-unstable configuration is half 
an order of magnitude less. The angle of attack is smaller in the configurations with 
relaxed static stability, and, as was the case for control law 1, A^ is also smaller in 
the configurations with relaxed static stability while A^^j b 

larger. ’ 

When comparing the turbulence responses of control law 2 with those of control 
law 1, it is observed that while A^^, A^^j -^Tj increase for control law 2, 

the remaining turbulence responses decrease. Because angle-of-attack and pitch-rate 
feedback are present for both sets of results and because the short-period damping is 
greater for all configurations and Mach numbers for control law 2, the smaller A val- 
ues are the result of the increase in short-period damping. The tail angle of attack 
increases even though the responses of some of its component parts decrease. The 
center-of-gravity vertical accelerations for control law 2 are greater than those for con- 
trol law 1 as a consequence of the corresponding decrease in pitch angle. The torque 
increased approximately 10 percent in each configuration over the torque for control 
law 1, but the bending moment decreased approximately 5 percent. Again, at both Mach 
numbers, the tail loads for the nominal -unstable configuration are significantly larger 
^at least 20 percent for A-p; almost 100 percent for than the tail-load responses 

for the baseline configuration, and the tail-load responses for the reduced-unstable con- 
figuration are significantly less ^60 percent reductions in A-p; 25 percent reductions in 
AM,b) than the tail-load responses for the baseline configuration. However, as before, 
the tail-load coefficient responses are significantly larger than the corresponding coeffi- 
cients in the baseline configuration. 
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Control Law 3 

Control law 3 uses forward-speed and pitch-rate feedback to satisfy the maneuver- 
ability requirement and at the same time maintain satisfactory flying qualities. 

Dynamic characteristics .- Table VTII contains the values of K^ and feedback 
gains for control law 3, the values of the natural frequencies and damping ratios of the 
phugoid and short-period modes, and the values of the apparent elevator angles per g. 
Table vni confirms that the maneuverability requirement (decreasing the apparent ele- 
vator angle per g of the RSS configurations 25 percent below the actual elevator angle 
per g of the baseline configuration) has been satisfied: at both Mach numbers the appar- 
ent elevator angles per g are within 2 percent of their required values. Table Vin also 
indicates that larger amounts of Ku and are required for the reduced-unstable 
configuration because its smaller tail results in less elevator effectiveness, and smaller 
amounts are required at the higher Mach number because aerodynamic forces are greater 
for a given elevator deflection. 

Figure 8 shows the phugoid and short-period roots (corresponding to the phugoid 
and short-period parameters of table VIII) plotted in the complex plane. The roots of the 
phugoid mode for each configuration are very highly damped and are well within the phu- 
goid FQC region. The roots of the short-period mode for each configuration are also 
well behaved, except that for the nominal-unstable configuration at Mach number 0.50 the 
short-period FQC are not satisfied. Because of the similarity of the frequencies and 
damping ratios of the nominal-unstable configurations at the two Mach numbers, it was 
elected to proceed with the analysis and to calculate the turbulence responses even though 
the flying-qualities criteria for one of the RSS configurations had not been satisfied. 

Figure 9 contains a typical root locus for changes in and Kq for the reduced- 
unstable configuration at Mach number 0.50. The solid lines for this figure represent 
increases in Ky (with Kg = 0.4e) and the dashed lines represent increases in Kg 
(with Ky = 0). When both gains are zero, as with the two previous control laws, the air- 
plane is statically unstable. As Kg increases, three things happen: (1) the dead-beat 

root at -0.024 in figure 9(a) combines with the servo dead-beat root (which was originally 
far to the left and out of the figure) to form the augmented short-period mode; (2) the roots 
of the third mode approach the real axis and split into two dead-beat roots; and (3) the 
divergent dead-beat root approaches, but does not enter, the left half plane. With K^ 
constant at 0.46, K^ is increased, as seen in figure 9(b), until the divergent dead-beat 
root enters the left half plane and combines with a root from the third mode to form the 
augmented phugoid mode. The short-period roots and the other root from the third mode 
remain unchanged with increases in Ky. In contrast to control laws 1 and 2 in which 
variation of K^^ only was sufficient to stabilize the configurations with RSS, control 
law 3 required changes in both Ku and Kg to stabilize such configurations. Also, 
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unlike control laws 1 and 2 (in which the location of the servo root was well removed from 
the locations of the phugoid and short -period roots) , the servo root for control law 3 is 
located "between" the phugoid and short-period roots. Such a servo root location, as 
stated in the section "Analysis," complicates the definition of flying-qualities criteria 
because of the resulting coupled motion. 

Turbulence responses .- Table IX contains the turbulence responses for all configu- 
rations and flight conditions for control law 3. When examining the turbulence responses 
from configuration to configuration at a given Mach number, it is observed that the 
forward-speed response increases significantly for the configurations with RSS. Angle 
of attack A^, for the configurations with RSS decreases on the order of 20 percent with 
respect to the baseline configuration primarily because of the increase in short-period 
damping. As with the tail loads calculated for the previous two control laws, the tail 
loads for the nominal-unstable configuration are larger than those for the baseline con- 
figuration, whereas the tail loads for the reduced-unstable configuration are smaller. 
However, the horizontal-tail loading intensity (as indicated by the tail-load coefficients) 
is larger than that for the baseline configuration for both RSS configurations and both 
Mach numbers. When comparing the turbulence responses of control law 3 with those of 
control law 2, it is observed that the elevator deflection is consistently smaller, and tail 
angle of attack, center-of-gravity acceleration, and torque are consistently larger. Bend- 
ing moment is larger for the nominal-unstable configurations but equal or smaller for the 
reduced-unstable configurations. 

Turbulence Response Comparisons 

Figure 10 contains bar charts of selected turbulence responses for each control law 
in terms of configuration changes. The responses presented are the center-of-gravity 
vertical acceleration, and torque and bending- moment coefficients. The letter designa- 
tions below each bar chart identify the configuration each response pertains to; B for 
baseline, NU for nominal unstable, and RU for reduced unstable. 

For all control laws and both flight conditions, the center-of-gravity vertical accel- 
eration and both tail-load coefficients are larger for both RSS configurations than for the 
baseline configuration. With only one exception, the center-of-gravity vertical accelera- 
tion and torque coefficient for the nominal-unstable configuration are consistently larger 
than similar responses for the reduced-unstable configuration for the same control laws 
and flight conditions. The exceptions occur for control law 3 at a Mach number of 0.50 
(fig. 10(c)). With only one exception the bending- moment coefficients for the reduced- 
unstable configuration are consistently larger than the bending -moment coefficients for 
the nominal -unstable configuration. The exception occurs for control law 3 at a Mach 
number of 0.75 (fig. 10(c)). In all cases the bending-moment coefficients for the reduced- 
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unstable configuration are more than twice as large as the bending-moment coefficients 
for the baseline configuration for the same control laws and flight conditions. 

Figure 11 contains bar charts of selected turbulence responses for both RSS con- 
figurations in terms of control law changes. The responses presented are center-of- 
gravity vertical acceleration, and torque and bending-moment coefficients. The number 
designations below each bar chart identify the control law that pertains to each response. 

For both configurations and both flight conditions, the center -of-gravity vertical 
acceleration and torque increase with increasing control law designation. The center-of- 
gravity vertical acceleration increases on the order of 20 percent and the torque increases 
on the order of 25 percent. This trend indicates that for the gains chosen, increasing 
maneuverability (control laws 2 and 3) is more costly, at least in terms of acceleration 
and torque, than matching the flying qualities of the baseline configuration (control law 1). 
For the nominal-unstable configuration (fig. 11(a)), the bending moment is also largest 
for control law 3 , but the responses are almost identical (within 5 percent) for control 
laws 1 and 2. For the reduced-unstable configuration (fig. 11(b)), however, the bending 
moment is largest for control law 1 and has almost identical responses (within 4 percent) 
for control laws 2 and 3. 


CONCLUDING REMARKS 

An analytical study has been conducted to determine airframe and horizontal-tail- 
load turbulence responses of control-configured rigid vehicles with relaxed static stability 
(RSS). The study included three vehicle configurations: one baseline configuration and 
two representative control-configured-vehicle configurations with RSS. The configurations 
with RSS relied on similar static stability augmentation systems (SSAS) for stability and 
each SSAS featured angle -of- attack (a) feedback, pitch-rate (0) feedback, and forward- 
speed (u) feedback to the elevator. The flight conditions chosen for the study were Mach 
numbers 0.50 and 0.75 at an altitude of 6100 m. Three control laws were examined which 
involved satisfying stability requirements, flying-qualities requirements, and maneuver- 
ability requirements by varying the SSAS feedback gains. One control law (using a- and 
0-feedback) matched the baseline short-period flying qualities exactly with no regard for 
maneuverability considerations. The other two control laws (one using a- and 0-feedback; 
the other usii^ u- and 0-feedback) increased the apparent maneuverability of the RSS con- 
figurations over that of the baseline configuration and at the same time satisfied the 
flying-qualities criteria. 

The results of this study indicated that 

1. Configurations with relaxed static stability can be stabilized either with angle- 
of-attack and pitch-rate feedback to the elevator or with forward-speed and pitch-rate 
feedback to the elevator. 
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2. For both flight conditions and each control law, the center-of-gravity vertical 
acceleration and the horizontal-tail-load intensity were larger for either RSS configura- 
tion than for the baseline configuration. 

3. For both RSS configurations, using the SSAS to increase maneuverability was 
more costly (on the order of 5 percent) in terms of center-of-gravity vertical accelera- 
tion and tail torque, but less costly (on the order of 5 percent) in terms of tail bending 
moment, than using the SSAS to match the baseline flying qualities. 

4. For both RSS configurations, using forward-speed and pitch-rate feedback to 
increase maneuverability was more costly (on the order of 10 percent) in terms of center- 
of-gravity vertical acceleration and tail-load intensity than using angle-of-attack and 
pitch -rate feedback to increase maneuverability. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
January 9, 1976 
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APPENDIX A 


FLYING-QUALITIES CRITERIA 

The flying-qualities criteria (FQC) described in reference 15 are defined in terms 
of upper and lower limits of the stability characteristics (damping ratios and undamped 
natural frequencies) of the phugoid and short-period modes. The purpose of this appen- 
dix is to illustrate an alternate way of presenting the FQC such that they are independent 
of Mach number. 

Damping-ratio and undamped-natural-frequency limits are conveniently represented 
in the complex plane. Constant damping ratio is represented by a pair of radial lines 
through the origin of the complex plane and symmetric with respect to the real axis; con- 
stant undamped natural frequency is represented by a circle whose center is the origin of 
the complex plane. Allowable limits of damping ratios and undamped natural frequencies 
are therefore represented by allowable areas (that is, FQC regions) in the complex plane 
which are bounded by radial lines and circular arcs. For convenience, the real and 
imaginary axes of the complex plane will have the units at* and icot*, respectively. 


Phugoid Flying Qualities 

Because the phugoid is such a long-period mode, its flying qualities are specified 
in terms of damping-ratio limits only. For the particular level, category, and classifi- 
cation of the present study, the damping-ratio limits are (ref. 15) 



C g 0.04 (Al) 

which corresponds to phugoid roots located 
within the (unshaded) FQC region. The 
phugoid FQC region in sketch (b) is the same 
for any Mach number. 
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Short-Period Flying Qualities 

Short-period flying qualities are specifi^ed in terms of both damping-ratio and 
undamped-natural-frequency limits. For the particular level, category, and classifica- 
tion of the present study, the dampii^-ratio limits are (ref. 15) 

2.00 s ^ > 0.30 (A2) 

The undamped-natural-frequency limits are functions of the quantity n^/a, which is 
expressed in reference 17 as 


0! K 1 

Mg-iZsMq 


(A3) 


After nondimensionalizing the aerodynamic coefficients and simplifying the resulting 
expression, equation (A3) becomes 


^z _ ^ 

a mg ^Lof 


(A4) 


For the particular level, category, and classification of the present study, reference 15 
presents a figure similar to the following slwtch (sketch (c)) for specifying the undamped- 
natural-frequency limits 



Sketch (c) 

The curves defining the upper and lower frequency limits are straight lines, each with a 
slope of +1/2 on the log-log plot. Such a slope suggests the following equations for cal- 
culating the upper and lower limits of ct^n,sp terms of n^/cx 
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(A5) 



(A6) 


where and C 2 are constants. When equation (A4) is substituted into equations (A5) 
and (A6), the upper and lower limits of Wn,sp expressed as 



= CiUo 


pS 

2mg 




(A7) 


(^n,sp)£ - Cl" (A8 

Equations (A7) and (A8) indicate that for an airplane of fixed geometry and fixed weight 
flying at constant altitude (which is the case in this study), a change in flight Mach num- 
ber will result in changes in Uq and Cl^. The changes in Uq and Cl^ result in 
different upper and lower limits of ct>n,sp each Mach number. 

Multiplying both sides of equations (A7) and (A8) by t* ^“ 2 ^^ expresses the fre- 
quency limits in terms of undamped natural reduced frequencies 


(kn,sp)^ - K,sp)^^t* - Cl 


(A9) 


(kn,sp)£ - (^n,sp)/* " ^2 f 

Equations (A9) and (AlO) indicate that for an airplane of fixed geometry and fixed weight 
flying at constant altitude, a change in flight Mach number will result in changes in Cj_^^ 

only. The upper and lower limits of gp at one Mach number may be expressed in 
terms of the upper and lower limits of k^^gp at another Mach number by 


' pS 
2mg 


■^a 


(AlO) 
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(All) 





(A12) 


In the present study the square roots of the ratio of the lift-curve slopes in equations (All) 
and (A12) are very close to unity (0.93), and consequently the upper and lower limits of 
kr. cr. at the two Mach numbers of this study are very nearly equal. Therefore, for sim- 
plicity and consistent with other assumptions made in the study, the upper and lower lim- 
its of kn sp th® two Mach numbers are assumed to be equal and resulted in the single 
short-period (unshaded) FQC region for both Mach numbers shown in sketch (d). 


io)t* 



at* 


Sketch (d) 


The short-period FQC region in sketch (d) is the same for any Mach number. 
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EFFECT OF FEEDBACK GAINS ON MANEUVERABILITY 

The purpose of this appendix is to define maneuverability and to show the effect of 
feedback gains and on the apparent elevator angle per g (the maneuverability 

parameter of this study). 


Maneuverability 

The cases of actual and apparent maneuverability for conventional (statically stable) 
airplanes and for airplanes with relaxed static stability are considered in this section of 
the appendix. Maneuverability in the longitudinal degrees of freedom usually is defined 
by the elevator angle per g, which is the elevator deflection required to trim the airplane 
in steady accelerated flight (such as a steady pull-up maneuver) of constant load factor. 
This deflection A6/An is related to the control-stick deflection per g (A6s/An) through 
the gear ratio Kg 



An 


(Bl) 


Sketch (e) illustrates the deflections of equation (Bl) for a statically stable configuration 
and sketch (f) illustrates the deflections for a configuration with relaxed static stability 
without a static stability augmentation system (SSAS). 



Sketch (e) 


Linkage 



Sketch (f) 
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For a configuration with relaxed static stability and with an SSAS, a part of the 
required elevator deflection is provided by the SSAS, which permits the control stick to 
be deflected in the conventional sense. This condition is described by the following 
equation: 


A 6 _ ^ ^^s 


+ K^ 


An ° An 


(B2) 


Equation (B2) is illustrated in sketch (g), where the pilot- commanded deflection 
ASj = Kg A 6 g and the SSAS- commanded deflection A 62 = Kq, Aa + A0 



Sketch (g) 


For both the conventional configuration without an SSAS, illustrated in sketch (e), 
and the configuration with relaxed static stability with an SSAS, illustrated in sketch (g), 
the control-stick deflection is in the conventional sense for a steady pull-up maneuver. 
From equation (Bl), the control-stick deflection per g is 

^ M (B3) 

An Kg An 

and is a measure of the actual maneuverability. From equation (B2), the control-stick 
deflection per g is 

^ = J_f^ . Kc, ^ - Ka^^ (B4) 

An Kg\An An ° Any 

where 

- K Ag y. Ag\ _ as’ 
yAn An ^ Any An 
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so that equation (B4) takes the form 

^^s 1 A5' 

An Kg An 


(B5) 


The quantity A6'/An is an apparent elevator deflection, which is taken to be a measure 
of apparent maneuverability. If feedback gains and are zero, equations (B4) 
and (B5) are measures of the actual maneuverability. 

The actual elevator angle per g, A6/An in equation (B3), is expressed as 


An 



(B6) 


Equation (B6) is obtained by algebraically manipulating a similar expression in refer- 
erence 12. The quantities Cl^q and are stability derivatives; f^-y^c is the 

nondimensional tail length; and hm, hn, and h are the locations of the maneuver point, 
neutral point, and center of gravity in terms of fractions of the wing mean aerodynamic 
chord. The quantities hj„ - h and hn - h are the maneuver margin and static margin, 
also in terms of fractions of the wing mean aerodynamic chord. 

The apparent elevator angle per g, A6'/An in equation (B5), may be expressed in 
terms of an apparent maneuver point and an apparent neutral point and expressed in a 
form similar to equation (B6) 


Ml 

An 



(B7) 


The quantities hj^ and h^ are the apparent maneuver point and the apparent neutral 
point and the quantities hjj^ - h and h^ - h are the apparent maneuver margin and the 
apparent static margin, all in terms of the wing mean aerodynamic chord. It remains to 
be shown that the apparent maneuver and neutral points are functions of the actual maneu- 
ver and neutral points and feedback gains and K^. 

Equations of Motion 

The equations of motion (eq. (1)), including the control law (eq. (2) and the last row 
in eq. (1)), are simplified and manipulated in this section of the appendix. Expressions 
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for the apparent neutral point, the apparent maneuver point, and, ultimately, the apparent 
elevator angle per g are derived from these simplified and manipulated equations later 
in the appendix. 

The equations are simplified by employing the short-period approximation (which 
eliminates the u variable) and by making the calm-atmosphere assumption (which elimi- 
nates Q!g). An additional simplification is made in the control law, which is written here 
in terms of circular frequency and with the u and 0!g terms removed 


6 = 


1 + iojtj^L 


K„(a) + K 


giicoB) 


(B8) 


It is assumed that at the frequencies of interest (the short period having the highest), the 
quantity icotj is sufficiently less than unity so that it may be set to zero. Written now 
in terms of the Laplace variable, the simplified control law is 


6 = K„(o!) + K^(se) 


(B9) 


Equation (B9) now becomes the last row of the equations of motion. As a consequence of 
the short-period and calm-atmosphere assumptions, the equations of motion have been 
reduced from a system of four simultaneous equations (in variables u, a, 6, and 6) 
with gust input on the right-hand side to a system of three simultaneous equations (in 
variables a, 6, and 6) with no gust input. 

This system of three simultaneous equations (in variables a, 6, and 6) is manip- 
ulated into a system of two simultaneous equations (in variables a and B) in the follow- 
l).ng manner: the simplified control law (eq. (B9)) is removed from the system of equations, 
Vhich leaves two equations in three variables; equation (B9) is then substituted for 6 in 
the remaining two equations and the following system of two simultaneous equations (in 
variables a and 6) results: 
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Apparent stability derivatives within equation (BIO) are defined in terms of the stability 
derivatives and feedback gains as 



(BU) 

= ^Zq + 

(B12) 


(B13) 

C!v, = Crv. + KA -4’ 

^mq '-mq . ♦ 

(B14) 


When equations (Bll) to (B14) are substituted into equation (BIO), it becomes 
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(B15) 


The equations of motion expressed by equation (B15) are now in proper form for deriving 
expressions for the apparent neutral point, the apparent maneuver point, and the apparent 
elevator angle per g in terms of feedback gains Kq, and K^. 

Apparent Neutral Point 

The derivation of the expression for hn as a function of and a discussion of 

the effect of on h^ are presented in this section. The actual neutral point appears 

within the following expression (ref. 12) for CniQij the static stability derivative. 


Cmo; 



(B16) 
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Similarly, the apparent neutral point appears within the expression for the 

apparent static stability derivative, 


C* 


a 



(B17) 


Substituting equation (B16) into equation (B13), equating equations (B13) and (B17), and 
then solving for h^ gives the following expression for the apparent neutral point in 
terms of the actual neutral point and feedback gain 

^mr 

hn = hn - 77-^ (B18) 

CLq, 

From equation (B18) it is clear that when is zero, the apparent and actual neutral 

points are the same. However, as increases, the apparent neutral point moves aft; 

this movement is stabilizing for a configuration which is initially statically unstable. 

Figure 12 contains plots of the apparent neutral point and the apparent static margin 
as functions of for both statically unstable configurations and both flight conditions. 

The originally statically unstable configurations become statically stable at values of 
for which the apparent static margin becomes positive. Because the Cnig stability 
derivative is directly proportional to the area of the horizontal tail (all other parameters 
constant), the slope of h„ as a function of Kq, is also directly proportional to the area 
of the horizontal tail (all other parameters constant). This dependence on tail area 
accounts for the shallower slopes on the plots of h^ against Kq, for the reduced- 
unstable configuration. It also accounts for the larger values of Kq, required for sta- 
bility for the reduced-unstable configuration. 

Apparent Maneuver Point 

The derivation of the expression for hj^j as a function of Kq, and K^ and a 
discussion of the effects of Kq, and K^ on hj[j^ are presented this section. 

From the expansion of the stability determinant of equation (B15) the following 
expression is satisfied only when the short-period mode is neutrally stable (that is, when 
the center of gravity is at the maneuver point) (ref. 12) 


^ZQ,<^mq - 


(B19) 


When equations (Bll), (B13), and (B14) are substituted into equation (B19), the following 
expression results: 
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prj ~ (B20) 

When equation (B16) and the following approximation (ref. 12) are substituted into equa- 
tion (B20), 


^Zq! = ~^L, 


a 


the following equation results: 


+ CmqCzgKa + CzgCm^KojK^ 


- 2/iCL^fh - hnj + 2/iCmgK, 




(B21) 


Equation (B21) is satisfied only when the center of gravity is located at the maneuver point 
hjjj or, in this case, the apparent maneuver point h^- Substituting hj^ for h in 
equation (B21) and solving for h^^ gives the following expression for the apparent 
maneuver point: 




— \ hn ” 


n 


2 /i j 2 /iCl 


’mqCzgKo, 


2MCmgKc, + C 




a 


Kq-Kq ^ 


(B22) 

The actual maneuver point is expressed as a function of the actual neutral point as 
(ref. 12) 


^m ~ hjj 



(B23) 


When equation (B23) is substituted into equation (B22), the apparent maneuver point is 
expressed as a function of the actual maneuver point and feedback gains K^, and 




m 




^mn^Zc^a “ 2/iCjn Kqj + C 2 -C 


a 




a 


CmgKg p 


(B24) 
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Equation (B24) indicates that when and are both zero, the apparent and actual 

maneuver points are the same. However, as or or both increase, the apparent 

maneuver point moves aft. 

Figure 13 contains plots of the apparent maneuver point and the apparent maneuver 
margin as functions of Kq, and for both statically unstable configurations and both 
flight conditions. The solid curves represent variations in with equal to zero 

and the dashed curves represent variations in with Kq, equal to zero. A surface 
would be necessary to show the effects of varying and simultaneously. Fig- 
ure 13 shows that the apparent maneuver point is more sensitive to than to and 

that the slopes of the curves are shallower for the reduced-unstable configurations and 
reflect the dependence of h^ on tall area. 


Apparent Elevator Angle per g 

A discussion of the effects of K^, and on A5’/An is presented in this 
section. 

The expression for the apparent elevator angle per g In terms of the apparent 
maneuver point and the apparent neutral point was presented in equation (B7) and 
repeated here 


as’ _ *^L,o ^m ~ ^ 


(B25) 


Equations (B18) and (B24) of the previous two sections show that the apparent neutral 
point h^ is a function of feedback gain and that the apparent maneuver point h^ 

is a function of feedback gains Kq, and Kq. Thus, the apparent elevator angle per g is 
a function of feedback gains and K^. Because h^ and h^ reduce to hm and 

hjj when both feedback gains are zero, the apparent elevator angle per g reduces to the 
actual elevator angle per g when both feedback gains are zero. 

Figure 14 contains plots of the apparent elevator angle per g as fxinctions of 
and for both statically unstable configurations and both flight conditions. The solid 
curves represent variations in Kq, with equal to zero and the dashed curves 
represent variations in with equal to zero. A surface would be necessary to 

show the effects of varying and simultaneously. Figure 14 shows that the 

apparent elevator angle per g is more sensitive to than to and that, contrary 

to figures 12 and 13, the slopes of the curves show very little dependence on tail area, but 
show considerable dependence on Mach number. A comparison of figures 13 and 14 shows 
that, as predicted, the apparent elevator angle per g becomes negative at the same values 
of gain at which the apparent maneuver margin becomes positive. 
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TABLE I.- PHYSICAL DIMENSIONS OF EXAMPLE CONFIGURATIONS 



Configuration 

Quantity 

Baseline 

Nominal 

unstable 

Reduced 

unstable 

Mass of airplane, m, 

5 730 

5 730 

5 444 

Pitch moment of inertia, Iyy> 

35 766 

34 953 

35 419 

Wing span, b, m 

13.54 

13.54 

13.54 

Wing mean aerodynamic chord, c, m 

2.55 

2.55 

2.55 

Wing area, S, m2 

31.8 

31.8 

31.8 

Sweep angle of wii^ quarter chord, A 25 , deg 

28.6 

28.6 

28.6 

Horizontal -tail span, bt, m 

5.35 

5.35 

3.78 

Horizontal-tail mean aerodynamic chord, ct, m 

1.47 

1.47 

1.04 

Horizontal-tail area, St, m2 

7.16 

7.16 

3.58 

Sweep angle of horizontal-tail 40 percent chord, At, deg . . . 

26.5 

26.5 

26.5 

Tail lei^th, £t> ni 

5.10 

4.85 

4.54 

Elevator chord ratio 

0.35 

0.35 

0.35 

Elevator servo characteristic time, tj, sec 

0.066 

0.066 

0.066 


TABLE n.- CENTER-OF-GRAVITY POSITION AND STABILITY 
PARAMETERS (IN FRACTION c) 

[Altitude = 6100 m] 


Quantity 

Mach number, 0.50 

Mach number, 0.75 

Baseline 

Nominal 

unstable 

Reduced 

unstable 

Baseline 

Nominal 

unstable 

Reduced 

unstable 

h 

0.32 

0.42 

0.31 

0.32 

0.42 

0.31 

hn 

.38 

.38 

.27 

.37 

.37 

00 

hn-h 

.06 

-.04 

-.04 

.05 

1 

b 

CJI 

-.03 

hm 

.40 

.39 

.28 

.39 

.39 

.29 

hm - h 

.08 

-.03 

CO 

0 

1 

.07 

CO 

0 

1 

-.02 
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TABLE in.- STABILITY DERIVATIVES OF EXAMPLE CONFIGURATIONS 


[Altitude = 6100 m] 


Stability 

derivative 

Mach number, 0.50 

Mach number, 

0.75 

Baseline 

Nominal 

unstable 

Reduced 

unstable 

Baseline 

Nominal 

unstable 

Reduced 

unstable 

Cx 

■^u 

-0.040 

-0.040 

-0.040 

-0.061 

-0.061 

-0.061 

Cxq, 

.102 

.102 

.100 

.038 

.038 

.037 

^L,o 

.216 

.216 

.205 

.096 

.096 

.091 


-.072 

-.072 

-.068 

-.124 

-.124 

-.117 

Cz 

-4.80 

-4.80 

-4.66 

-5.54 

-5.54 

-5.39 

Cz . 

-1.65 

-1.57 

-.733 

-2.22 

-2.11 

-.989 


-3.46 

-3.29 

-1.54 

-4.02 

-3.81 

-1.78 


-.502 

-.502 

-.251 

-.583 

-.583 

-.291 

Cmu 

0 

0 

0 

0 

0 

0 

Cm 

IHq, 

-.287 

.191 

.186 

-.276 

.276 

.161 

Cm- . 

mQ, 

-3.29 

-2.97 

-1.30 

-4.44 

-4.01 

-1.76 

Cmq 

-6.91 

-6.23 

-2.73 

-8.02 

-7.24 

-3.17 

Cmg 

-1.00 

-.953 

-.446 

-1.16 

-1.11 

-.517 
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TABLE IV.- STABILITY AND MANEUVERABILITY PARAMETERS FOR CONTROL LAW 1 

[Altitude = 6100 m] 


Mach 

number 

Configuration 

Feedback gains 

Short period 

Phugoid 

ASVAn, 

deg/g 

rad/rad 

^ 9 ’ 

rad/rad/sec 

ko 




0.50 

Baseline 

0 

0 

1.73 X 10-2 

0.560 

6.80 X 10-4 

0.058 

-0.968 

.50 

Nominal unstable 

.44 

.03 

1.70 

.557 

6.26 

.062 

-.952 

.50 

Reduced unstable 

.88 

.15 

1.73 

.557 

5.99 

.070 

-1.896 

.75 

Baseline 

0 

0 

1.61 

.653 

3.39 

.198 

-.333 

.75 

Nominal unstable 

.41 

.03 

1.61 

.640 

2.77 

.242 

-.327 

.75 

Reduced unstable 

.61 

.11 

1.63 

.652 

2.63 

.267 

-.634 


TABLE V.- TURBULENCE RESPONSES FOR CONTROL LAW 1 
[Altitude = 6100 mj 




Feedback 

gains 

Turbulence responses 

Mach 

number 

Configuration 

Ka, 

Kg, 


Aa, 






Atj«, 

AM,b> 


^M,b’ 



rad 

rad 

1 

rad 

rad 

rad 

rad/sec 

rad 

g units 

N-m 

N-m 

1 

1 



rad 

rad/sec 

m/sec 

m/sec 

m/sec 

m/sec 

m/sec 

m/sec 

ra/sec 

m/sec 

m/sec 

m/sec 

m/sec 

0.50 

Baseline 

0 

0 

3.77 X 10"3 

6.04 X 10-3 

1.45 X 10-3 

6.20 X 10-3 

3.45 X 10-3 

0 

6.76 X 10-2 

19.6 

81.4 

2.28 X 10-4 

9.45 X 10-4 

.50 

Nominal unstable 

.44 

.03 

2.97 

6.00 

1.60 

4.69 

3.28 

1.32 X 10-3 

7.45 

22.2 

143.7 

2.58 

1.67 X 10-3 

.50 

Reduced unstable 

.88 

.15 

3.23 

6.04 

1.50 

4.99 

3.12 

2.46 

7.28 

7.34 

67.6 

2.41 

2.23 

.75 

Baseline 

0 

0 

6.56 X 10-4 

3.97 

7.71 X 10-4 

2.10 

2.59 

0 

1.12 X 10-1 

27.7 

109.9 

1.43 

5.68 X 10-4 

.75 

Nominal unstable 

.41 

.03 

3.23 

3.90 

9.19 

1.16 

2,77 

8.17 X 10-4 

1.29 1 

33.8 

217.1 

1.75 

1.12 X 10-3 

.75 

Reduced unstable 

.61 

.11 

4.49 

3.90 

8.23 

1.45 

2.13 

1.09 X 10-3 

1.22 

10.5 

85.4 

1.53 

1.25 



TABLE VI.- STABILITY AND MANEUVERABILITY PARAMETERS FOR CONTROL LAW 2 


[Altitude = 6100 m] 


Mach 

number 

Configuration 

Feedback gains 

Short period 

Phugotd 

A6'/An, 

deg/g 

rad/rad 

Kg, 

rad/rad/sec 





0.50 

Baseline 

0 

0 

1.73 X 10-2 

0.560 

6.80 X 10-4 

0.058 

-0.968 

.50 

Nominal unstable 

.25 

.11 


1.029 

3.18 

.161 

-.725 

.50 

Reduced unstable 

.60 

.03 

9.98 X 10"2 

.619 

5.99 

.057 

-.718 

.75 

Baseline 

0 

0 

1.61 X 10-2 

.653 

3.39 

.198 

-.333 

.75 

Nominal unstable 

.27 

.06 


1.024 

8.91 X 10-5 

.630 

-.253 

.75 

Reduced unstable 

.45 

.02 

9.47 X 10-3 

.682 

2.80 X 10-4 

.244 

-.253 


TABLE vn.- TURBULENCE RESPONSES FOR CONTROL LAW 2 
|Altitude = 6100 nij 




Feedback 

gains 

Turbulence responses 

Mach 

number 

Configuration 

^ a > 


aq. 

Aq ,, 


>1 

A^, 

Ag, 

A^n’ 

Atp, 

AM,b> 

Ac^> 




rad 

rad 

1 

rad 

rad 

rad 

rad/sec 

rad 

g units 

N-m 

N-m 

1 

1 



rad 

rad/sec 

m/sec 

m/sec 

m/sec 

m/sec 

m/sec 

m/sec 

m/sec 

m/sec 

m/sec 

m/sec 

m/sec 

0.50 

Baseline 

0 

0 

3.77 X 10-3 

6.04 X 10-3 

1.45 X 10-3 

6.20 X 10-3 

3.45 X 10-3 

0 

6.76 X 10-2 

19.6 

81.4 

2.28 X 10-4 

9.45 X 10-4 

.50 

Nominal unstable 

.25 

.11 

4.99 X 10'^ 

5.51 

1.73 

8.14 X 10-4 

1.68 

9.19 X 10-^ 

7.74 

23.6 

136.1 

2.75 

1.59 X 10-3 

.50 

Reduced unstable 

.60 

.03 

2.59 X 10“3 

5.91 

1.64 

3.84 X 10-3 

1.94 

1.85 X 10-3 

7.68 

7.87 

61.8 

2.58 

2.03 

.75 

Baseline 

0 

0 

6.56 X 10"^ 

3.97 

7.71 X 10-4 

2.10 

2.59 

0 

1.12 X 10-^ 

27.7 

109.9 

1.43 

5.68 X 10-4 

.75 

Nominal unstable 

.27 

.06 

3.51 

3.67 

9,91 

1.10 

2.30 

6.59 X 10-4 

1.34 

36.0 

213.1 

1.86 

1.10 X 10-3 

.75 

Reduced unstable 

.45 

.02 

2.89 

3.87 

8.89 

1.02 

1.63 

8.63 

1.27 

11.2 

80.5 

1.64 

1.18 
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TABLE Vm.- STABILITY AND MANEUVERABILITY PARAMETERS FOR CONTROL LAW 3 


(Mtitude = 6100 m] 


Mach 

number 

Configuration 

Feedback gains 

Short period 

Phi^oid 

A6'/An, 

deg/g 

Ku, 

rad/m/sec 

Kg, 

rad/rad/sec 

ko 



c 

0.50 

Baseline 

0 

0 

1.73 X 10-2 

0.560 

6.80 X 10-4 



.50 

Nominal unstable 

- 1.3 X 10-4 

.30 

4.79 

.823 

7.82 X 10-5 


-.718 

.50 

Reduced unstable 

- 2.7 

.46 

2.26 

.948 

6.28 

.833 

-.738 

.75 

Baseline 

0 

0 

1.61 

.653 

3.39 X 10-4 

.198 

-.333 

.75 

Nominal unstable 

- 6.2 X 10-5 

.18 

4.94 

.719 

5.94 

.797 

-.254 

.75 

Reduced unstable 

- 7.5 

.20 

2.38 

.896 

4.99 

.858 

-.244 


TABLE IX.- TURBULENCE RESPONSES FOR CONTROL LAW 3 
[Altitude = 6100 mj 


Feedback gains Turbulence responses 


nSSber Configuration 














Ku, 

K^, 

Au, 

A^, 

Aof^, 

Ag, 

Ag, 

A§, 

Aarj 

Aq^, 

^M,b’ 


Ap , 

^M,b’ 


rad 

rad 

1 

rad 

rad 

rad 

rad/sec 

rad 

g units 

N-m 

N-m 

1 

1 




m/sec 

rad/sec 

m/sec 

m/sec 

m/sec 

m/sec 

m/sec 

m/sec 

m/sec 

m/sec m/sec 

m/sec 

m/sec 

0.50 

Baseline 

0 

0 

3.77 X 10’3 

6.04 X 10'3 

1.45 X 10-3 

6.20 X 10-3 

3.45 X 10-3 

0 

6.76 X 10-2 

19.6 

81.4 

2.28 X XO'4 

9.45 X 10-4 

.50 

Nominal unstable 

- 1.3 X 10-4 

.30 

4.10 

4.99 

1.97 

4.76 

2.63 

7.41 X 10-4 

8.53 

26.7 

143.7 

3.10 

1.67 X 10-3 

.50 

Reduced unstable 

- 2,7 

.46 

6,17 

4.69 

2.00 

7.35 

3.06 

1.36 X 10-3 

8.89 

9.43 

61.8 

3.10 

2.03 

.75 

Baseline 

0 

0 

6,56 X 10-4 

3.97 

7.71 X 10"4 

2.10 

2.59 

0 

1.12 X 10-1 

27.7 

109.9 

1.43 

5.68 X 10-4 

.75 

Nominal unstable 

- 6,2 X 10-5 

.18 

1.56 X 10-3 

1 

3.35 

1.14 X 10-3 

4.17 

3.44 

5.77 X 10-4 

; 1.48 

40.8 

226.0 

2.11 

1.17 X 10-3 

.75 

Reduced unstable 

- 7.5 

.20 

1.84 

3.25 

1.07 

4.82 

3.16 

6.10 

1.43 

13.2 

80.1 

1.93 

1.17 
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(a) Short-period region. (b) Phi^oid region. 


Figure 2.- Flyii^ qualities criteria regions. 
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(a) Short-period roots. 
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(b) Phugoid roots. 

Figure 4.- Roots of longitudinal modes for control law 1. 
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(a) Short-period locus. 
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Third mode 



(b) Phugoid locus. 

Figure 5.- Typical root locus for control law 1. Reduced-unstable 

configuration; M = 0.50. 
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(a) Short-period roots. 
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(b) Phugoid roots. 

Figure 6.- Roots of longitudinal modes for control law 2. 
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(a) Short-period locus. 
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(b) Phugoid locus. 


Figure 7.- Typical root locus for control law 2, Reduced-unstable 

configuration; M = 0.50. 
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(a) Short-period roots. 
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(b) Phugoid roots. 

Figure 8.- Roots of longitudinal modes for control law 3 
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(a) Short-period locus. 
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(b) Phugoid locus. 

Figure 9.- Typical root locus for control law 3. Reduced-unstable 

configuration; M = 0.50. 
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Figure 12o- Effect of angle-of-attack feedback gain on effective neutral point 
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Figure 14.- Effects of angle-of-attack feedback gain and pitch-rate feedback 
gain on effective elevator angle per g. 
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proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Defails on the avaUabilify of these pubiications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


Washington, D.C. 20546 




